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Abstract. The gravitational microlensing experiments in the direction of Large Magellanic 
Cloud (LMC) predict a large amount of white dwarfs (~ 20%) filling the galactic halo. How- 
ever, the predicted white dwarfs have not been observed at the galactic halo. To interpret the 
microlensing results and resolving the mentioned problems, we use the hypothesis of spatially 
varying mass function of MACHOs, proposed by |Kerins &: Evans (1998)] (hereafter KE). How- 
ever the KE model is not compatible with the duration distribution of events [(Rahvar 2004a) | 
Here we use more realistic power-law model of MF, dn/dm oc m° for the MACHOs of halo. 
The index of MF in this model changes from —2.7 for stars with m > IMq at the central part 
of galactic halo to the substellar regime with an upper limit of —1 at the edge of halo. We 
show that in contrast to the abundant brown dwarfs of galactic halo, heavy MACHOs can be 
responsible for the microlensing events in the direction of LMC. 



1. Introduction 

The rotation curves of spiral galaxies and the Milky Way show that these type of 
galaxies have dark halo component. The most trivial candidate for the dark halo structure 
is the baryonic matter that can be in the form of gaseous or MAssive Compact Halo 
Objects (MACHOs). Since MACHOs are expected to be too light to be luminous and 
difficult detectable, [Paczyhski (1986) proposed an indirect method so-called gravitational 



microlensing to observe them indirectly. Following his suggestion several experiments such 
as EROS and MACHO started monitoring millions of LMC stars for one decade and 
observed less than 20 events in this direction fAlcoc k et al. 200?11 ILasseree et al. 2000|l . 
Due to the degeneracy nature of gravitational microlensing problem, it is impossible to 
obtain the mass, distance and transverse velocity of the lenses by measuring the duration 
of events. The only way is using statistical studies, comparing the microlensing events 
with the models. The outcome of this study is the mean mass of MACHOs and their mass 
contribution in the galactic halo. For the standard Galactic model with the Dirac-Delta 
mass function (MF), halo is comprised by 20% of MACHOs with the mass of 0.5Mq. 
Comparing the luminous mass of Milky Way with that of halo, we can conclude that the 
white dwarfs of halo should have twice mass of the ordinary stars of the disk and bulge. 
KE proposed using Dirac Delta spatially varying MF as a solution to this problem but it 
has been shown that the duration distribution in this model is not compatible with the 
LMC microlensing candidates. Rahvar (2004a) used a likelihood analysis to find the best 
parameters for the model. Here we use more realistic MF for the MACHOs of halo as 
the power law function where the mass scale and the index of MF changes monotonically 
from the center to the edge of galaxy. We show that duration distribution of events in this 
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model is compatible with the LMC microlensing candidates. The effective mean mass of 
lenses, as the mean mass of observed lenses in this model is larger than the mean mass of 
overall lenses of the galactic halo. This means that the galactic halo with the dominant 
brown dwarfs can produce microlensing events compatible with the observed data. 
The organization of paper is as follows: In Section 2 we introduce the hypothesis of 
spatially varying MF of KE and power-law models. In Section 3 we compare the LMC 
microlensing events of MACHO experiment with the expected distribution from the 
model and section 4 contains the conclusion of this work. 



2. Spatially Varying MF 

The tradition in the microlensing analysis is using Dirac-Delta MF for the MACHOs 
of galactic halo. This function has been used for the sake of simplicity to apply it to 
an arbitrary MF of halo. However this approach hides the information about the contri- 
bution of the MACHOs with various masses in the gravitational microlensing candidates. 
Star and baryonic cluster formation theories (e.g., Ashman 1990 Carr 1994; Dc Paol is et al. 1995ll 
predict the variation of MF in interstellar medium depending on the density of clouds 
forming the stars. Shadmehri (2004) and Elmegreen (2004) also showed that instead of 



universal initial mass function we can have a composite with various contributions from 
the brown dwarfs, solar masses and high mass stars, depending on the density of en- 
vironment. Since the halo density changes as ~ l/'r^, we can conclude that the MF of 
MACHOs may change by distance from the center of galaxy. The inner halo comprises 
partly visible stars, in association with the globular cluster population, while the outer 
halo comprises mostly low-mass stars and brown dwarfs. 

The hypothesis of using spatially varying MF has been proposed by KE. They used a 
spatially varying MF as 

MF{r) = 5[M - M{r)], (2.1) 
where the mass scale in this model decreases monotonically as a function of r as 

Mc/(^)'-/«--, (2.2) 

Mjj and Ml are the upper and lower limits for the mass of MACHOs and Rhaio is the 
size of galactic halo. In the CMD plus baryonic halo models, Rhaio the size of halo which 
contains MACHOs can be extended more than the halo size. 

Here we use more realistic spatially varying MF than the KE model. The MF of stars 
is taken as the power low function dn/dm oc m", where the index of function depends 
on the mass scale of stars. The upper limit of index (Chabrier 2002) for the substellar 
regime can be —1 and the lower limit (Scalo 1986) of —2.7 for the stars massive than the 



solar mass. 

Here we use the following anzats for the spatially varying MF: 

dm Tnolr) ' 



(2.3) 



where m(r) > nio(r) and the index of power-law changes as a(r) = 1-1-^ 2.7. 

The typical mass of MACHOs, mo(r) changes as M[/(-^)''/^''"'° . Here we take the 
Mu = IMq, Ml = O.OOIMq and Rhaio = lOOfcpc. Figure n shows the distribution of 
mass of MACHOs as a function of distance from the center of galaxy. 
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Figure 1. The mass distribution of MACHOs as a function of a; = Dd/Dos, wiiere Doi and Dos 
are tlie observer-lens and observer-source distances. The upper panel and lower panels show the 
distributions of mo and mass of MACHOs as function of distance. 

3. Microlensing Analysis in the Spatially Varying MF 

Here in this section our aim is to compare the microlensing data with the theoretical 



model. First we calculate the theoretical and experimental optical depths IjAlcock et al. 2000|l 
and compare them with each others. Using the experimental efficiency of observation, 
the observed optical depth obtain as follows: 



Tobl 



AE 



(3.1) 



where ti is the duration of ith event and E is the exposure time of observation which 
for the 5.7 yrs observation of LMC by MACHO experiment, E is 6.12 x 10^. The optical 
depth results from 13 LMC microlensing candidates is about 4.43 x 10~*. One the other 
hand, using the power-law MF for the MACHOs of galactic halo in the standard halo 
model, we can produce the rate of microlensing events ^ . The theoretical optical depth 
is calculated as 



T, 



expected 



= - j -emt, 



(3.2) 



where we consider hundred percent of halo is made by MACHOs and e{t) is the obser- 
vational efficiency. The ratio of the observed to the theoretical optical depth results the 
MACHOs contribution in the galactic halo, / = — — . Here in this model 23 percent 

'^expected 

of halo is made by MACHOs. Figure |2] shows the distribution of duration of events in the 
standard halo model, considering the observational efhciency of MACHO experiment. 
We do a Monte-Carlo simulation to show quantitatively how the contribution of massive 
lenses in the microlensing events is more than the light ones. According to the distribu- 
tion of matter along the line of sight the lenses are selected by the probability function 
as follows: 

dP 

— oc p{x)y^x{l - x), (3.3) 

where p is the matter density and x is the lens-source to the observer-source distances. 
After selecting the position of lenses, we use the equation (|2.3|) to obtain their masses. The 
duration of events are generated through te = and compared with the observational 
efficiency to be selected or rejected. Re is the Einstein radius of a lens and vt is its 
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Figure 2. The duration distribution of LMC microlensing events {te) , considering the observa- 
tional efficiency of MACHO experiment in Standard halo model and spatially varying power-law 
MF. 

transverse velocity with respect to our line of sight. The mass of selected events are used 
to obtain the mean mass of lenses. We called the mean mass of these lenses as the active 
mean mass. In contrast to the active mean mass we define passive mean mass of lenses 
as the mean mass of overall lenses in the halo. The active and passive mean masses of 
lenses in this model are obtained as the 0.40 and 0.05 solar masses. This means that 
brown dwarfs are the abundant MACHOs of halo while white dwarfs are responsible to 
the microlensing events. 

4. Comparison with the observed data 

In this section we compare the expected duration distribution from the model with 
duration distribution of the observed LMC candidates. Two statistical parameters as the 
mean and the width of the duration distribution of events IjGreen and Jedamzik 20021 
IRahvar 2004b|l are used for our comparison. The width of the duration distribution for 
the Nobs'th observed candidate is defined by 

AtE = Max[j=i^N,^^]{tj) - AIin[j^i^ff^^^]{tj). 

The mean and the width of duration distribution of MACHO candidates are 97 and 188.2 
days. Figure 13 compares the observed mean and width of duration of events with what 
we expect form the theoretical model. We perform a Monte-Carlo simulation to generate 
the mean and the width of distribution of events. Thirteen events is picked up from the 
duration distribution of events (Figure^) and each time the mean and the width of events 
are calculated in each set to produce their distributions. Comparing these distributions 
with the observation show that the standard halo model with the power-law spatially 
varying MF is in good agreement with the LMC data. 

5. Conclusion 

In this work we use an anzats for the spatially varying power-law MF for the MACHOs 
of galactic halo. The advantage of using this MF is that (i) it is more realistic than Dirac- 
Delta MF, (ii) the expected duration distribution of events is in good agreement with 
the LMC microlensing data and (iii) the active mean mass of lenses as the mean mass of 
lenses contributing in the microlensing events obtained about 0.40Mq compared to the 
passive mean mass of lenses as the mean mass of overall lenses of the galactic halo with 
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Figure 3. Upper and lower panels show the theoretical distributions of the mean and the 
width of microlensing events in standard halo model with spatially varying power-law MF. The 
cross signs are the mean and the width of duration of observed LMC candidates of MACHO 
experiment. 

0.05Mq . We showed that in a halo with abundant brown dwarfs, massive MACHOs with 
the white dwarf sizes can be responsible for the microlensing events that we are observing 
in the experiment. 
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